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ABSTRACT 

Detailed knowledge of the mesopheric sodium layer characteritics is crucial to estimate 
and optimize the performance of Laser Guide Star (LGS) assisted Adaptive Optics 
(AO) systems. In this paper, we present an analysis of two sets of data on the meso- 
spheric sodium layer. The first set comes from a laser experiment that was carried out 
at Cerro Tololo to monitor the abundance and altitude of the mesospheric sodium in 
2001, during six runs covering a period of one year. This data is used to derive the 
mesospheric sodium column density, the sodium layer thickness and the temporal be- 
havior of the sodium layer mean altitude. The second set of data was gathered during 
the first year of the Gemini MCAO System (GeMS) commissioning and operations. 
GeMS uses five LGS to measure and compensate for atmospheric distortions. Analy- 
sis of the LGS wavefront sensor data provides information about the sodium photon 
return and the spot elongation seen by the WFS. All these parameters show large 
variations on a yearly, nightly and hourly basis, affecting the LGS brightness, shape 
and mean altitude. The sodium photon return varies by a factor of three to four over 
a year, and can change by a factor of two over a night. In addition, the comparison of 
the photon returns obtained in 2001 with those measured a decade later using GeMS 
shows a significant diff'erence in laser format efficiencies. We find that the temporal 
power spectrum of the sodium me an altitude follow s a lin ear trend, in good agreement 
with the results reported by Pfro mmer fc Hicksonl (|201(l ). 
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1 INTRODUCTION 

1.1 Mesospheric sodium layer 

Sodium (Na) atoms are believed to be deposited in the 
high atmosphere by meteoritic ablation. They form a layer 
whose average altitude lies in the 90-95 km altitude range 
above sea level. This altitude can vary from 85 km to 
105 km depending on the particular physical conditions. 
The sodium abundance displays typical seasonal variations 
of a factor of 2 to 4, with the minimum and maximum 
seasonal a bundance occurring in summer and winter, re- 
spectively (|Moussaoui et al.l I2OI0I ). Significant abundance 
variations on hourly, daily and yearly time scales have 
been repo rted, eve n on ti me scales of a few seconds to a 
few hours llClemesh a & Ta kahashi|[l99 5: Milon ni et al.ll 19981: 
O'Sullivan et all I2OOQ : .Butler et al.l I200Q : .Michaille et all 
200 ll ). The sodium layer mean altitude and width also ex- 
hibit significant variations on short and long time scales: the 



layer is in average lower in altitude (by 1 to 2 km) and thin- 
ner (by 1 km) at the e quinoxes, and high er in altitude and 
thicker at the solstices (|Papen et al.lll99a ). Most of the short 
time scale abundance and mean altitude variations can be 
traced down to the appearance of relatively short-lived, high 
density, thinner layers with in the main sodium la yer, which 
are refered to as sporadics (|Clemesha et alJI 19961 ). 



1.2 Laser Guide Stars in Adaptive Optics 

Observations with Adaptive Optics (AO) are limited to cer- 
tain areas on the sky due to the requirement of a stellar 
source to measure the wavefront distortions. A solution to 
this sky cover age problem is to creat e an artificial guide star 
with a laser (|Fov fc Labevriel Il985l ). A so-called "Sodium 
LGS" can be created from resonant backscattering of meso- 
spheric sodium atoms. Today, almost a ll leading telescope s 
are equipped with LGS-AO systems (|Amico et al.ll2010l l. 
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and recently, the Gemini Observatory produced the first 
sodium LGS constellation to feed a Multi-Conjugate Adap- 
tive Optics (MCAO) system at the Gemini South 8-meter 
telescope. With the increasing prevalence of such LGS sys- 
tems, it is becoming more critical to understand and investi- 
gate the intrinsic properties and characteristics of the meso- 
spheric sodium layer. For instance, a good knowledge of the 
sodium abundance is crucial to understand the fiuctuations 
of the brightness of the LGS during the year, and hence the 
impact and limitations on the AO performance. This infor- 
mation is also helpful to optimize queue observations. Also 
of importance is the study of the sodium layer interaction 
with a given laser format so this interaction can be opti- 
mized, as different laser formats can lead to very di fferent 
resuhs (jRochester et al.ll2012l : iHolzlohner et al.ll2010l l. 

1.3 Sodium parameters for Laser Guide Stars 

The following characteristics of the sodium layer are key for 
the AO performance: 

• Mesospheric sodium column density. 

• Temporal behavior of the sodium layer mean altitude. 

• Sodium layer thickness. 

Mesosphenc sodium column density 

The mesospheric sodium column density is the most crucial 
parameter to characterise so that the effectiveness of the 
sodium LGS technique can be maximised. This is because, 
in absence of saturation or at low saturation levels typically 
achieved with currently existing LGS AO systems, the 
required laser power per beacon is directly proportional to 
the sodium column density. Therefore, the measurement of 
the minimum sodium column density over a year sets an 
upper limit for laser power requirements. Additionally, the 
seasonal and nightly statistics of sodium density fluctua- 
tions, which are equivalent to LGS magnitude fluctuations, 
provide insights into the expected long and short term 
system performance such that the efficiency and flexibility 
of queue scheduling can be maximised. 

Temporal behavior of the sodtum layer mean altitude 
The rate of mean sodium altitude variations is of prime 
importance for optimum removal of the focus mode in AO- 
corrected science images. Focus adjustments in the LGS 
path of AO systems must distinguish atmosphere-induced 
focus terms from slow drifts of the guide star altitude. The 
faster the mean sodium altitude varies, the more difficult it 
is to distinguish those two effects and adequately correct for 
them. GeMS is equipped with a Slow Focus Sensor (SFS) 
that monitors the sodium altitude drifts. Knowledge of the 
required rate at which this sensor should run is thus very 
important. On the other hand, for AO systems working with 
a LGS constellation like GeMS, the differential focus error 
between the LGS is also of importance. A differential focus 
between the LGS would produce a signal that could not 
be properly treated by a tomographic reconstructor, and 
would therefore lead to non-negligeable loss of performance 
out of the AO system. The temporal behavior of this 
differential error is key for the desig n of future AO s y stems 



on Extremely Large Telescopes i Pfrommer et al.1 |2012| : 
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iHerriot et aLll2012l ; iDiolaiti et al]|2012l ). 



Sodium layer thickness 

The sodium thickness will determine the LGS spot elonga- 
tion as viewed from a subaperture located near the edge of 
the telescope pupil (this is ~1 arcsec for an 8 m telescope 
and comparable to the LGS spot size). The larger the tele- 
scope, the more elongated the spot, the more sensitive the 
adaptive optics system is to noise and therefore the more 
laser power is required to achieve a given performance. 



1.4 Sodium data 

With the increasing use of sodium LGS in astronomical AO, 
many studies have been done in the past years to character- 
ize more precisely the properties of the sodium layer in the 
light of AO requirements. Models have been developed to 
predict the sodiu m abundance and expec t ed photon return 
for d ifferent sites (jMoussaoui et al" ] |2009l : IHolzlohner et al.l 
I2OIOI 1. The temporal behavior of sodium profiles, and the 
variatio ns of the mean altit u de have also been studied i n 
detail (|Herriot et al.1 l2006l : iPfrommer fc HicksonI I2OI0I '). 
However, whereas the sodium layer had been observed 
quite extensively at several locations, no measurements had 
been made yet in Chile at or near the latitude of Cerro 
Pachon. The Gemini Observatory therefore initiated a 
year-long sodium monitoring campaign at the Cerro Tololo 
Inter-American Observatory (CTIO) located only a few 
kilometers away from the Gemini South telescope. In 2001, 
during a series of 6 runs, a laser experiment was carried 
out at Cerro Tololo to monitor the abundance and altitude 
of the mesospheric sodium. The goal was to characterize 
yearly, nightly and hourly the m esospheric sodium varia - 
tions at CTIO/Gemini latitude (Id'Orgeville et al.1 bOOSl ). 
This set of data is compared with the results obtained 
during the first year of the Gemini MCAO System (called 
GeMS) commissioning. GeMS started on-sky commission- 
ing in January 2011 , at a rate of 5 to 7 nights per month 
fRig aut et all l2012i ). This first period of commissioning 
lasted 5 months, after which GeMS entered a shutdown 
phase for engineering upgrades. On-sky operations resumed 
in November 2011, and continued up to May 2012. During 
this first period of commissioning, we gathered data about 
the sodium return as seen by the LGS wavefront sensors, 
and the differential focus between the LGS. The differential 
focus error analysis will be presented in a forthcoming paper. 

This paper is organized as follows. Section 2 presents 
the CTIO campaign along with the data reduction and re- 
sults. Section 3 presents the GeMS data, data reduction and 
results. Section 4 compares the results of both data sets, and 
discusses the impact of the laser format on performance. 



2 THE CTIO/TOLOLO CAMPAIGN 
2.1 Overview 

The experimental se tup has been described in detail in 
Id'Orgeville et al.l (|200al here we only present its main char- 
acteristics. The experiment involved launching a laser beam 
whose wavelength was tuned to the sodium D2 absorption 
line (589 nm) to the sky such that mesospheric sodium 
atoms were excited to higher energy levels. Sodium atoms 
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Table 1. Data summary, "useable data" means that the pho- 
tometry was reliable enough for laser profiles not to depend on 
weather variations (e.g. variable cirrus cover). 



Run number 


Date 


Night with useable 
laser data 


1 


Feb, 11-20, 2001 


7 


2 


May, 2-11, 2001 


3 


3 


Aug 31 - Sep 6, 2001 


2 


4 


Nov 25 - Dec 1, 2001 


6 


5 


Feb 23 - Mar 2, 2002 


5 



were excited with a low-power continuous-wave laser whose 
interaction with so dium atoms was relatively well-known 
l|Milonni et al.lll999l ). The laser equipment included a 6-7 W 
multi-line argon-ion laser pumping a commercial ring-dye 
laser. The laser output power was monitored in real time at 
a 0.3 Hz rate so that sodium density fluctuations could be 
calibrated out from laser power fluctuations. The on-sky 589 
nm launched power was in the ~ 100-200 mW range and sent 
to the sky by a custom launch telescope. The laser beam, a 
truncated gaussian of ~250 mm at the 1/e^ intensity points, 
was projected about 3 degrees off zenith at a fixed angle. 
The 0.9 m CTIO telescope and the University of Michigan 
Schmidt telescope were used to image the laser guide star 
thus created. They were respectively at a distance of 140 
m and 110 m away from the laser projection system yield- 
ing a spot elongation of ~35 arcsec. Images were taken with 
successive exposures of typically 10 seconds but this period 
of data acquisition varied depending on the photometry of 
the night. It reached up to 30 seconds on some of the obser- 
vations. Standard stars were also observed a few times per 
night in order to derive atmospheric transmission. Provided 
that all due telescope calibrations were performed, and that 
the atmosphere optical transparency was measured during 
the night, the CCD data then contained all information nec- 
essary to retrieve the sodium column density, sodium layer 
width and relative altitude. The absolute altitude could also 
be derived by a triangulation method based on the natural 
guide star trails seen on the CCD behind the laser streak, 
however this method has not been implemented on the data 
yet, hence only relative altitudes will be presented here. 

2.2 Data and data reduction 

Observations were spread over one year from February 2001 
to February 2002 to allow characterization of hourly, nightly 
and yearly variations of the sodium layer parameters. Five 
runs of 7 to 10 nights were performed in February, May, 
September, November 2001 and February 2002. Those 
particular months were chosen to match the expected 
minimum (November/December) and maximum (May) of 
the sodium column density sinusoid-like variations. Table [T] 
indicates the dates of each run and the resulting number of 
useful nights when data was taken. 

Data reductio n follows the procedure/method de- 
scribed in Id'Orgeville et all (|2003l '). The first part of the 
work consisted of retrieving, cleaning and organizing all the 
data from the compact disks (CDs) and notebooks the ten 
year-old data was stored on. All images were fiat-fielded 
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Figure 1. Example of reduced image obtained on the 0.9m tele- 
scope CCD, showing the laser streak. This image was taken at 
UT=00h51min on September 4th, 2001. 



and bias-substracted using regular IRAF procedures. The 
nature of the fiat fields varied from dome to sky fiats 
depending on the run and photometry of the night. In 
one instance (Run 4) dark currents were also substracted. 
Fig. [T] shows an example of the laser streak (data from May 
10th 2001) once the image has been reduced. The next 
steps were to calibrate the fiux received from the LGS into 
photons/s/cm'^/W and convert it into a sodium column 
density in atoms/cm^. A first part of this calibration 
process consisted in deriving a zero point (ZP) and the 
atmospheric transmission for each night. This was based 
on observations of standard stars conducted during the 
night. These standard stars were observed either at the 0.9 
m directly, or the nearby Schmidt telescope. Additionally, 
some standard stars were observed through the same Na 
filter, whereas other stars were observed with a V filter. 
We derived the different factors to convert the flux into a 
flnal ZP based on the notes and data that we had. The 
second part of this calibration concerns normalization by 
the laser fiux. Laser power data for the last three runs 
(September 2001, November 2001 and February 2002) had 
been logged automatically which minimised any transcript 
errors. However, laser power data for the first two runs 
(February 2001 and May 2001) was only available as 
handwritten notes, therefore increasing data uncertainty 
and limiting the calibration accuracy. The conversion from 
sodium return measured in photons int o sodium column 
density follows the same assumptions as in Id'Orgeville et al.l 
( 2003 ). using an absorp tion cross section value of 1.0 10^^ 
cm^ (see equation (1) in IdHDrgeville et al.l (|2003[ l and figure 
2 in lMilonni et al.l (| 19981 ) 1. Note that the elapsed time of 10 
years between data taking and data reduction unfortunately 
results in some uncertainties regarding calibrations. This 
is especially true for Run 1 and 2 during which the data 
log was only available as handwritten notes. The errors 
introduced by the uncertainties in the calibration data are 
hard to quantify due to the different observation conditions 
at different points of the year. It is our estimate that 
such errors, in a worst case scenario, should be no larger 
than 25%. We also expect the seasonal sodium abundance 
variations to have larger amplitudes than what can be 
accounted for by this error. Therefore, the error introduced 
by the uncertainties in the calibration procedure only adds 
to a noise floor that the sodium abundance variation will 
be clearly above, ensuring our proceeding analysis is sound. 
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Figure 2. Two sodium layer intensity profiles extracted from 
the same night; February 25th 2002, UT=01h35min (solid line) 
and UT=07h45 (dashed line). Flux is expressed in Arbitrary Unit 
(AU). 
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2.3 Results 

2.3.1 Sodium profiles 

By integrating over the image it is possible to derive an 
intensity profile that can be used to investigate such intrin- 
sic properties of the sodium laser as equivalent width and 
mean altitude. Fig. [2] shows two profiles that have been ex- 
tracted from the 25th of February 2002 (run 5) at respec- 
tively UT=01h35min (solid line) and UT=07h45 (dashed 
line). The former profile demonstrates a typical laser profile 
shape, while the latter profile shows that a double asymmet- 
ric peak has developed, where the intensity of the second 
peak is almost twice stronger as the first. This example il- 
lustrates the variability that one can expect during a given 
night and emphasizes the importance of characterizing such 
behavior. 

Concatenating all the different profiles, one can follow this 
evolution over a whole night. Profiles have been generated 
for all the nights which had data that was of high quality 
(see Table 1). Fig.[3]shows three different sample nights. The 
previous example of February 25th is shown as the last plot 
at the bottom on Fig. [S] The event seen on the 2nd profile 
above appears around 05h UT and it develops up to the end 
of the night. The different examples in Fig. [3] illustrate the 
high variabil ity of the sodium layer i n terms of profiles, as 
report ed by jPfrommer et al.l l|2010l ) ; jPfrommer fc HicksonI 
1I2OIOI) . 



2.3.2 Sodium equivalent width 

For each profile, we estimate an equivalent width for the 
sodium thickness. For this, we first normalize each profile 
by its maximum. The equivalent width is then found by 
forming a rectangle with a height equal to one, and finding 
the width such that the area of the rectangle is equal to the 
area in the profile. For profiles with irregular shapes, the 
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Figure 3. Sodium profiles for 3 sample nights. From top to bot- 
tom, the nights of: 4th September 2001, 27th November 2001, 
25th February 2002. The white columns correspond to time in- 
tervals in the observation night when no data was taken because 
there was a problem with the telescope and/or laser. 



equivalent width is more robust than a Full Width at Half 
Max (FWHM) estimation. The equivalent width gives an 
estimate of the spot elongation to be expected on the Laser 
Guide Star WaveFront Sensor (LGSWFS). In the small angle 
approximation, the angular elongation 7 is given by: 

' H^ + HAH ' ^ ' 

where AH is the thickness of the sodium layer, H is the 
low altitude of the sodium layer, and L is the distance be- 
tween the pupil subaperture and the laser launch telescope. 
The top plot of Fig. |4] shows the evolution of the equivalent 
width for the night of 25th February 2002. The bottom plot 
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Figure 4. Variation of the sodium profile equivalent width. Top 
is for the night of the 25th February 2002. Bottom is the aver- 
age equivalent width per night for all the data-set. Error bars 
represent the minimum and maximum values per night. 



of Fig.Ushows the average width per night for the entire set 
covering one year of data. Average equivalent widths are on 
the order of 10 km, as expected, which leads to a spot elonga- 
tion of ~1 arcsec for the Gemini telescope whith its on-axis 
launch configuration. Significant variations on hourly and 
daily basis are seen. The equivalent width can change by al- 
most a factor of two within a given night. A small correlation 
with seasons can be detected: the sodium layer is somewhat 
thinner (resp. thicker) around September (resp. February) 
by ~1 k m. This was a l so ob served at other latitudes as re- 
ported bv lPapen et al.l (| 19961 ). The largest equivalent widths 
seen on this data are on the order of ~15 km, which corre- 
sponds to an elongation of 1.3 arcsec for the Gemini tele- 
scope. To avoid clipping the LGS spots on the WFS sub- 
aperture, this maximal expected elongation should be used 
to set the requirements on the WFS subaperture Field Of 
View (FoV). The GeMS subaperture FoV is 2.8 arcsec, which 
allows for the spot elongation and the kernel due to seeing. 



2.3.3 Sodium mean altitude 

As described in the introduction, variations in the mean alti- 
tude of the sodium layer will be interpreted by the LGS AO 
system as focus variations, deteriorating the performance of 
the AO correction. A variation of Aft in the mean altitude 
of the sodium layer will cause an rms phase error (piston- 
removed) of: 



1 D^sin{ip) 
16\/3 {h-hoY 



Ah 



(2) 



where D is the diameter of the telescope, is the zenith 
angle, h is the mean sodium altitude and ho is the telescope 
altitude. Using Marechal's approximation, we can estimate 
the loss of performance in term of Strehl Ratio (SR) due to 
this defocus phase error as: 



SR = exp 



(3) 



where A is the observation wavelength. 



In the top plot of Fig. O we show an example of the 
derived mean altitude for a profile acquired during the 
night of the 28th of November 2011. For this given night, 
the mean altitude changed by ~6 km, which, according to 
eq. [2] and |3] would correspond to a loss of SR of ~100% 
in _ff-band (1.65 /im). The bottom plot shows how the 
defocus error changes over different periods of time At, for 
four sample nights. This plot is computed as an average of 
the mean altitude differences for a given At. The altitude 
difference is then converted into a defocus error by using 
Eq. [2] As the figure shows, the larger the time between 
re-focussing, the greater the variation in the the mean 
altitude, and thus the more the focus error grows. Note that 
the dashed lines represent a loss of 1%, 5% and 20% SR in 
-ff-band (1.65 nm). Note also the large variation from night 
to night: one can expect a factor of two in the amplitude 
of the defocus error. To cope with this issue, LGS-AO 
systems are usually equipped with an independent focus 
sensor looking at a natural guide star so as to detect any 
deviations from the optimal focus, since focus variations 
on a NGS can only be induced by the atmosphere. The 
rate at which this sensor must work sets the requirement 
on the guide star magnitude, hence on the sky coverage. 
For the examples of Fig. (5] if we do not want to degrade 
the SR by more than 5%, a focus update should be done 
every ~3 min. If we reduce the requirement to 1% {Ah 
= 43 m), the focus update should be done every ~30 s. 
In the current scheme for GeMS, we are using refocusing 
rates ranging from 1 second to 5 min depending on the 
magnitude of the NGS. Good focus measurements are 
typically obtained for Is exposure time on NGS with R 
< 13.0. This means that if we want to keep the focus 
error down to a loss of 1% of SR in 7J-band, our current 
limiting magnitude is about R = 16.7. For fainter stars, the 
focus error will grow. The current limiting magnitude is 
above the original requirement of R — 18.5, and thus this 
limitation in magnitude is impacting sky cove rage. This is a 
known issue of the current hardware in GeMS (|Neichel et al.l 
|2010D . Performance is expected to be improved in the future. 
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Figure 5. Top: Variation of the mean altitude of the sodium 
layer for the night of 28th of November 2011. Bottom: Average of 
the mean altitude change for given At in minutes. Four nights are 
illustrated here, from top to bottom: 28th November 2001, 27th 
November 2001, 4th September 2001 and 25th February 2002. 
The horizontal dotted lines represent a loss of 1%, 5% and 20% 
SR in //-band (1.65 fim). 



To generalize the results presented in Fig. [5] we com- 
pute the Power Spectral Density (PSD) of the mean altitude 
variations for all the nights in our sample. This i s shown in 
Fig. [6] As found bv lPfrommer fc HicksonI (|2010l ). the PSD 
of the mean altitude variations is well fitted by a model de- 
fined by: 



(4) 



where, for our set of data, we find a best fit by using a=35 
m'^Hz"^ and /?=-1.9 (dashed line in F ig. O. This is in good 
agreem ent with the results found by IPfrommer fc HicksonI 
l|2010l ) from data acquired in the northern hemisphere, 
where they derive a=30±20 m^Hz"^ and /3=-1.95 ± 0.12. 
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Figure 6. Power Spectral Density of the mean profile altitude. 
Dots are for individual data points, the dot-dash line is the aver- 
age over the different nights and the dashed line is the best fit to 
the data. 



2.3.4 Sodium return flux 

This parameter is critical for LGS AO systems as it is di- 
rectly related to the flux received by the LGSWFS of the AO 
system. The sodium return (or brightness of the LGS) can 
change rapidly during a night, and it is expected to reach its 
lowest average level during summer time, possibly in Novem- 
ber/December due to variations in the sodium abundance. 
In Fig. [71 we plot the evolution of the sodium photon return 
during the night of the 25th of February 2002 and for the 
night of the 28th of November. One can see that the return 
can vary by up to a factor of three within the same night. 
This impacts the AO correction, and should be taken into ac- 
coun t by automatic optimization procedures jNeichel et al.l 
l2010h . The seasonal variation observed during the one year 
duration of these measurements is shown in Fig [8] in terms 
of photons/cm^/s per Watt of projected laser power, given 
at the sodium layer (atmospheric transmission is not fac- 
tored in on the way back to the ground) and in terms of 
sodium column density. The average sodium column den- 
sity varies by a factor of ~2 between a summer minimum of 
3 10^ atoms/cm^ in November/February and a winter max- 
imum of 6 10® atoms/cm'^ in May. These results are in rea- 
sonable agreement with sodium abundance values reported 
in previous studi es (iPapen et al.lll996l : iMichaiUe et"aLll200ll : 
[m oussaoui et al. ||20lrf) and confirm that the sodium prop- 
erties for the Gemini-South site are similar to those of other 
sites. The impact of the seasonal variations on GeMS per- 
formance is discussed in section R3. 3. II 



3 GEMS COMMISSIONING 
3.1 Overview 

GeMS is the fi rst multi-sodium based LGS AO system used 
for a stronomv l|d'Orgeville et all |2008| . |2012| : iRigaut et all 
l2012h . The GeMS laser is a CW mode-locked 50 W laser. 
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Figure 8. Average photon return and Sodium abundance (in x 
10^ atoms/cm^) per night for all the nights that had useable data. 
The return is given at the sodium layer (atmospheric transmission 
is not factored in on the way back to the ground). Error bars 
represent the minimum and maximum values per night. 



and the projected beams are not controlled in polarization 
yet. Each LGS is seen by a dedicated LGSWFS made of a 
16x16 subaperture Shack-Hartmann lenslet array. Each sub- 
aperture is sampled by 2x2 pixels (quadcell configuration). 
The effective laser power projected into the sky is calibrated 
with respect to the GeMS laser system output power and the 
output power is monitored at a rate of 1 Hz. The projected 
beams are gaussian, with beam diameters on the order of 25 
cm at the 1/e^ intensity points. Spot size measured on the 
sky are ranging between 1.2 and 1.7 arcsec. LGSWFS data 
has been regularly saved during the GeMS commissioning 
nights. The photon return is extracted from this data. 



Figure 7. Variation of the sodium photon return for two nights: 
top is February 25th 2002 and bottom is November 28th 2001. 
The return is given at the sodium layer (atmospheric transmission 
is not factored in on the way back to the ground). 



whose spectro-temporal and spatial format is quite different 
from the Cerro Tololo laser format, and significantly less ef- 
ficient in exciting sodium atoms at equal power levels. The 
Gemini South (GS) laser contains two infra-red (IR) laser 
lines (respectively 1064 nm and 1319 nm), each created by 
one oscillator and multiple amplifiers, combined together in 
a lithium triborate (LBO) non-linear crystal. Each NIR laser 
is actively mode-locked, resulting in a 77 MHz pulse train 
with nominal pulse widths on the order of 300-400 ps. The 
spectral bandwidth of the 589nm laser has been measured 
to be on the order of 1.5-2.1 GHz. The laser beam is relayed 
from the output of the Laser system to the input of the 
Laser Launch Telescope (LLT) located behind the telescope 
secondary mirror, by a set of mirrors called Beam Transfer 
Optics (BTO). Five LGSs are produced after splitting the 
50 W laser beam into five lOW beams in the BTO. The 
throughput of the BTO has been measured to be ~50%, 



3.2 Data and data reduction 

GeMS started on-sky commissioning in January 2011, at a 
rate of 5 to 7 nights per month. Data was obtained during 
the first 5 months of 2011, after which GeMS entered a 
shutdown phase for engineering upgrades. On-sky opera- 
tions resumed in November 2011, and continued up to May 
2012 again at a rate of one run per month. We do not use 
the data from January 2011 and February 2011, because at 
that time the laser spots were not yet optimized, and their 
FWHM was larger than the LGSWFS field stop, so fiux 
measurements may be biased. From March on, laser spots 
were of the order of 1.3 to 1.5 arcsec, which is smaller than 
the field of view of the LGSWFS subapertures of 2.8 arcsec. 
For all these data points, the laser stabilization loop which 
is keeping each of the LGS in front of the LGSWFS was 
closed, so we do not expect fiux loses because of coupling 
with the subaperture FoV. Table [2] summarizes the data 
available for each run. All the data was taken at zenith, 
in the same conditions, and nights with clouds have been 
discarded. 

The data reduction is fairly straightforward in that 
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Table 2. Summary of data used for the sodium return character- 
ization from GeMS commissioning and average return computed 
for each of these data sets, as seen by the LGSWFS. 



Date 


Number of nights 


Sodium return 
(in ph/s/cm^/W) 


March 2011 


3 nights 


5.4 


April 2011 


4 nights 


8.7 


May 2011 


1 night 


9.5 


November 2011 


3 nights 


4.8 


December 2011 


5 nights 


3.5 


January 2012 


4 nights 


4.1 


February 2012 


4nights 


8.2 


March 2012 


5 nights 


6.5 


April 2012 


3 nights 


11.5 


May 2012 


2 nights 


13 



case. The photon return is measured at the LGSWFS level, 
by integrating the flux over all the subapertures. The flux 
is first converted into photons/second/cm^ by using the 
LGSWFS CCD detector gain. Then, the normalization by 
the laser power projected to the sky is applied. All the data 
is time-stamped, which facilitates its cross-correlation. The 
LGSWFS data is acquired at a frame rate ranging from 100 
Hz to 800 Hz. However, as the laser power is monitored 
at a rate of 1 Hz, we average the LGSWFS data over bins 
of 1 s. Finally, it is important to note that as the flux is 
measured at the LGSWFS level, it includes all the losses 
due to transmission of the atmosphere, the telescope itself, 
and of Canopus (the GeMS AO bench). To translate this 
flux in terms of photons/cm^ /s per Watt of projected laser 
power at the ground i.e. at the primary mirror of the receiv- 
ing telescope, one can use the following transmission: ~ 
0.8 (measured Gemini South telescope throughput), T^q ~ 
0.28 (measured AO bench and LGSWFS throughput at 589 
nm), and LGSWFS QE = 0.8 (detector quantum efficiency 
from constructor data). 



0) 
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3.3 Results 

3. 3. 1 Sodium return flux 

FigE] shows the sodium return measured by the LGSWFS 
for two sample nights: the 14th of November 2011, and the 
13th of December 2011. For these two nights, we see that 
the flux suddenly increases by more than a factor of two, 
which is most likely due to the presence of sporadics. This 
is something that is regularly seen, at a rate of once every 2 
or 3 nights, and it illustrates the high variability of the LGS 
brightness that can be expected. 

Fig. HQ] shows the return flux for all the data available. 
The amplitude variations over a night are represented by 
the errors bars that show the minimal and maximal return 
measured during each night. As presented in Fig. O we 
retrieve here the variations due to the sodium season. Lower 
return is seen in December, while higher return is around 
May. The difference between low and high season is on the 
order of a fact or of three to four , in good agreement wit h 
other studies (|Papen et al.l Il996l : iMoussaoui et al.l I2OI0I ). 
This information can be used to optimize the scheduling 
of laser operations: between May and September laser 
operations would be facilitated as far as photon return is 



Figure 9. Sodium return measured by the LGSWFS for 2 sample 
nights: 14th of November 2011 (top) and 13th of December 2011 
(bottom). 

concerned. Unfortunately, this seasonal variation of sodium 
abundance is correlated with the bad weather conditions 
and median seeing variations across the year. This means 
that when the sodium density is at its maximum, seeing 
is also at its maximum. This is limiting the use of AO 
instruments during the Chilean winter and for instance it 
is planned for GeMS to be shutdown during the months of 
July and August every year. Such large seasonal variations 
also directly impact GeMS performance, as the frame rate 
of the AO loop should be adjusted consequently in order 
to keep the level of photons measured by the LGSWFS 
more or less constant. In order to get a reasonable signal to 
noise ratio on the WFS, we have established that around 
35 ph/frame/pixels are needed. If running at 800 Hz, this 
translates into 10 ph/s/cm^/W at the LGSWFS. Looking 
at the results of Fig. 1101 we see that this condition is 
only reached during the higher sodium season, i.e. around 
May. During the months of November and December, we 
have to reduce the AO frame rate, sometimes as low as 



Characterization of the sodium layer at Cerro Pachon 9 




Table 3. Summary of the main characteristics for the Tololo and 
GeMS laser. 



200 300 
Time [days] 

Figure 10. Sodium return measured by the LGSWFS for all 
the data available over more than 1 year of commissioning. Error 
bars represent the minimum and maximum values per night. 



100 Hz in order to keep the right SNR at the LGSWFS level. 



4 GEMS AND TOLOLO SODIUM PHOTON 
RETURN COMPARISON 

This section presents a side-by-side comparison of the 
sodium photon returns measured in 2001-2002 above Cerro 
Tololo using a low power, 100 mW-class CW monomode 
laser (about 350 mW projected on sky), and in 2011-2012 
above Cerro Pachon using a high power, 50 W-class CW 
mode- locked laser (about 4 W projected on-sky per LGS). 
Cerro Tololo and Cerro Pachon are about 10.5 km away, 
close enough that the sodium layer properties can be as- 
sumed to be the same above these two sites, albeit for mea- 
surements made 10 years apart. Table [3] presents a summary 
of the Tololo and Pachon laser power and spectro-temporal 
formats, which emphasizes the fundamental differences be- 
tween the two laser probes. It is by now well-known that 
not all 589 nm lasers are created equal in their interac- 
tions with the sodium layer, and that key site character- 
istics such as the site latitude and the geo-magnetic field 
strength and its orientation also strongly influence sodium 
photon ret urn results for a g i ven on - sky laser pointing above 



that site ijRo chester et al 
Milonni et al. 1998: Michaille et al 



20121: 



Holzlohner et alj 

"I2OOI; Moussao ui et all 



ng above 

IT I2010I: 



2010l : ld'Orge^ille et al.li2003 : iMilonni et al.iil999 '). With this 



knowledge, we will compare relative sodium abundance vari- 
ations measured with both lasers, and discuss the relative 
efficiency of the Gemini South laser format with respect to 
the Tololo laser format. 

Fig. [H] presents sodium photon return results obtained 
on Tololo both in terms of photons/cm^ /s/W of projected 
laser power, at the sodium layer, and in terms of sodium col- 
umn density. Fig. [TO] presents sodium photon return results 
obtained on Pachon only in terms of photons/cm^ /s/W of 
projected laser power, at the GeMS LGSWFS. At this time 



Description 


Tololo Laser 


GeMS Laser 


Average Power 


150-350 mW 


35-45 W 


Type 


CW 


Pulsed 


Pulse repetition 




77 MHz 


Pulse length 




300-400 ps 


Spectral bandwidth 


0.5 MHz 


1.5-2.1 GHz 


Polarization 


Linear 


Elliptical (not controlled) 



there is no direct sodium abundance measurement concur- 
rent to our data to enable accurate calibration of the sodium 
column density corresponding to the GeMS sodium photon 
return data. This is unfortunate since it would have made 
possible to compare not only relative but also absolute vari- 
ations in sodium abundance. As it is, the conversion factor 
between sodium photon return from the GeMS laser and the 
corresponding sodium column density can only be estimated 
using the relatively rough assumptions presented below. 

In order to compare the Tololo and GeMS sodium pho- 
ton returns, we choose to express both of them in terms 
of photons/cm^/s per Watt of projected laser power at the 
ground i.e. at the primary mirror of the receiving telescope. 
For the GeMS data, we use the conversion factors intro- 
duced in Sect. 13.21 For the Tololo data, we use the atmo- 
spheric transmission derived in Sect. 12.21 Nightly sodium 
photon return averages obtained in 2001-2002 at Tololo (Fig. 
[S]) thus translate into equivalent returns of ~60 to ~130 
photons/cm^ /s/W at the ground, while nightly sodium pho- 
ton return averages measured with GeMS in 2011-2012 (Fig. 
nop are in the 13 to 70 photons/cm'^/s/W range at the 
ground. 

Fig. [S] shows that the monthly November/December 
2001 average for sodium abundance is on the order of ~3.5 
lO'' atoms/cm^, corresponding to ~65 photons/cm^ /s/W 
at the ground. Combining the monthly averages provided 
in Table [5] for the months of November and December 
2011 yields a weighted November/December 2011 average 
of ~10.5 photons/cm'^ /s/W at the ground. Assuming that 
the low season, November/December monthly average value 
in both data sets corresponds to the same low value for 
sodium column density monthly averages, this yields a con- 
version factor of about 3 between the Tololo CW laser for- 
mat and the GeMS laser format. Doing the same thing for 
high season values in May 2001 (Fig. [HI and May 2011/2012 
(Table[2} respectively yields ~116 photons/cm^ /s/W on the 
ground at Tololo and ~66 photons/cm'^ /s/W on the ground 
at Pachon, resulting in a smaller conversion factor of ~1.8. 
The high season approach is deemed to be much less reli- 
able based on what is known of typical sodium abundance 
variation behavior, where sodium abundance can spike at 
times due to high rates of sporadics over a set of nights, 
whereas low sodium column density values typically remain 
of the same order of magnitude. The real conversion fac- 
tor between the two laser formats therefore lies somewhere 
between 1.8 and 3, and is believed to be closer to 3. Not 
surprisingly, the Tololo monomode CW laser appears to be 
several times more efficient in exciting sodium atoms than 
the Gemini South laser by a large fraction. How large this 
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factor really is remains to be determined with higher accu- 
racy. 

The sodium photon return results for the Tololo laser 
can a lso be compared to table 2 presented in (| Telle et alj 
|2006| ). Telle's Fasor laser is a high power (50W) CW 
monomode laser for which results extrapolated to low 
("zero") power basically apply to the Tololo laser format. 
Telle reports that, based on extrapolating the Fasor on- 
sky measurements, such a laser format provides 150 (resp. 
140, 185, 115) photons/cmVs/W "at zero power" of pro- 
jected Fasor power at the ground at the Starfire Optical 
Range in New Mexico, USA, on October 12 (resp. Octo- 
ber 14, November 16, and December 22), 2005, for an es- 
timated sodium column density of 7.5 (resp. 7.0, 8.4, and 
5.9) 10^ atoms/cm^. Scaling these values down to our as- 
sumed, Tololo low sodium column density value of 3.5 10^ 
atoms/cm^ yields sodium photon return values of 70 (resp. 
70, 77, and 68) photons/cm^ /s/W of projected Tololo laser 
power at the ground, in rather good agreement with our es- 
timated value of ~65 photons/cm^/s/W. This comparison 
gives us reasonable confidence in the Tololo results presented 
in section 2 and in particular in the sodium column density 
values we have derived from them. 

Sodium abundance measurements have been performed 
on Cerro Pachon by the neighboring Andes LIDAR Observa- 
tory (AL O, the Universit y of Illinois sodium LIDAR experi- 
ment, see lSwenson et al] |200^)) during various observation 
runs since 2009. Unfortunately no ALO data is available con- 
current with GeMS sodium photon return data which would 
have permitted a direct calibration of the sodium abundance 
based on GeMS results. Some data is available in between 
runs though, which we can use to determine the validity of 
our earlier estimates of how much more efficient the Tololo 
laser is with respect to th e GeMS laser. ALO sodium abun- 
dance data (Liu A]|2012| 1 was obtained over the nights of 
January 28-February 2, 2012 (average sodium column den- 
sity ~3.9 10^ atoms/cm^) and March 20-23, 2012 (average 
sodium column density ~4.9 10^ atoms/cm^), in between 
the GeMS January 7-12, February 10-13, and March 10-14, 
2012 nights used to derive the sodium photon return results 
presented in Table 1. Table |4] presents a summary of the 
relevant ALO, Tololo and GeMS data which can be used to 
compare actual vs. extrapolated GeMS sodium photon re- 
turns in the January-March 2012 period. Expected sodium 
photon returns are calculated for the Tololo and GeMS laser 
formats assuming a conversion factor of 3 between these two 
formats. Measured and extrapolated GeMS results in the 
January-March 2012 period are qualitatively consistent, re- 
inforcing our belief that a factor of 3 is probably close to the 
mark. 

We are hopeful that, in the future, closer coordination 
with the Andes LIDAR Observatory will permit concur- 
rent GeMS laser and LIDAR laser propagation above Cerro 
Pachon. Direct calibration of the sodium abundance at the 
same time as GeMS sodium photon return data is gathered 
will likely prompt a review or at least a refinement of the 
conclusions presented here. 



5 CONCLUSION 

We have presented results on the sodium photon return 
and characteristics of LGS created above Cerro Pachon and 
Cerro Tololo based on two sets of data. Both were gathered 
over more than a year of observations and have been used 
to compute the sodium photon return per Watt of projected 
laser, the sodium profiles and the variations of the sodium 
layer mean altitude. The results show that: 

• A LGS exhibits a large variety of fiuctuations in shape 
and brightness. These changes can appear on very short (few 
seconds) or seasonal time scales. 

• Average sodium equivalent width are on the order of 
~10 km, which leads to a spot elongation of ~1 arcsec for 
an 8 m telescope using an on-axis launch configuration. 

• A small correlation of the sodium equivalent width with 
seasons is detected: the sodium layer is thinner by 1km 
around September. 

• The PSD of the sodium mean altitude variation is well 
fitted by a model defined by PSD(!/) = av^ . We derived 
Q=35 m^Hz"^ and /3=-1.9 in good agreement with previous 
data published in Northern hemisphere. 

• A refocusing rate of <~30 s is required for GeMS. 

• The sodium photon return shows a large variability over 
a year, with seasonal fiuctuations on the order of 3 to 4. 

• Nightly fiuctuations of photon return by a factor of two 
are regularly seen. 

• Comparing sodium photon returns obtained with two 
different lasers clearly illustrates the impact of the laser for- 
mat on coupling efficiency with the sodium layer. For the 
two laser format described in Tab. |3] and used in this study, 
we measure a difference in efficiency by a factor of ~3. 
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